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Abstract 

Bursting is a periodic transition between a quiescent state and a state of repetitive 
spiking. The phenomenon is ubiquitous in a variety of neurophysical systems. We 
numerically study the dynamical properties of a normal form of subcritical Hopf os- 
cillator (at the bifurcation point) subjected to a nonlinear feedback. This dynamical 
system shows an inhnite-period or a saddle-node on a limit cycle (SNLC) bifurcation 
for certain strengths of the nonlinear feedback. When the feedback is time delayed, 
the bifurcation scenario changes and the limit cycle terminates through a homo- 
clinic or a saddle separatrix loop (SSL) bifurcation. This system when close to the 
bifurcation point exhibits various types of bursting phenomenon when subjected to 
a slow periodic external stimulus of an appropriate strength. The time delay in the 
feedback enhances the spiking rate i.e. reduces the interspike interval in a burst and 
also increases the width or the duration of a burst. 
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1 Introduction 



When a system's activity alternates between a quiescent state and a state 
of repetitive spiking, the system is said to exhibit bursting behaviour [1,2,3]. 
These large excursions exist because the system is close to bifurcations from 
rest to oscillatory states. There are two bifurcations which characterize a burst; 
one which leads a rest state to repetitive spiking and the second which leads 
a spiking state to the rest state. As there are varieties of such bifurcations 
one finds a number of different types of bursts [1]. The bursting is ubiquitous 

1 E-mail: gautam@ipr.res.in 

2 E-mail: abhijit@ipr.res.in 



Preprint submitted to Elsevier Science 



7 February 2008 



in a variety of physical and biological systems, especially in neural systems 
where it plays an important role in communication among neurons [4,5,6,7]. 
The spiking rate within a burst and also its duration are considered important 
for the reliability of neuronal communication. 

In this paper we describe a mathematical model which contains the necessary 
ingredients to not only produce different types of bursts but to also control 
the width and the spiking rate within a burst. This becomes possible due to 
introduction of a time delay in the self-feedback term of the system. We may 
point out here that the self-feedback can mock up a variety of physical effects. 
In a collection of coupled systems, the term can be regarded as a source term 
representing the collective feedback due to the rest of the systems. A similar 
single equation model as a paradigm for studying the collective dynamics 
(synchronization, death states etc.) of an ensemble of coupled oscillators has 
been proposed and discussed in depth in [8]. The time delays, arising from 
finite conduction/propagation speeds, can certainly be significant in situations 
where the systems communicate over some distance. The feedback is time 
delayed to account for such propagation delays of the signals. 



2 The Model 

A topological normal form of the subcritical Hopf oscillator is given as: 



where z = x + iy is a complex variable and \x is the bifurcation parameter. 
The frequency of the oscillations is determined by u and b\z\ 2 . The parameter 
b which is called the shear parameter, determines how the frequency depends 
on the amplitude of the oscillations. Canonical models derived from the above 
and similar normal forms have been widely used in the past for studying 
the bursting behaviour of neurons [9,10,11]. In what follows, we alter the 
above model to incorporate a nonlinear feedback term and then study its 
repercussions on the bifurcation dynamics and consequently on the excitability 
of the system. 

The system dynamics near Hopf bifurcation point is of intertest so for our 
present studies we assume the oscillator to be residing at the subcritical Hopf 
bifurcation point i.e. = 0. We add a nonlinear feedback term to the equation. 
The feedback is time delayed to incorporate the unavoidable time delays in 
signal transmission in any physical or biological system. In our model we have 
introduced a quadratic nonlinear feedback which is the simplest and the lowest 
order nonlinearity that provides an excitable behaviour to the model as we 
see later. We denote the feedback strength by k and the time delay by r. The 
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modified oscillator equation can then be written as: 



z(t) = (i(u + b\z{t)\ 2 ) + \z(t)\ 2 - \z{t)\ 4 )z(t) - kz 2 {t - r) 



(2) 



We solve and carry out bifurcation analysis of this model equation rewritten in 
x and y coordinates. The next section gives the different bifurcation scenario 
with the nonlinear feedback strength k and the time delay r. 
Most of the simulations in the present paper have been carried out using soft- 
ware package XPPAUT [12] which also provides an interface to the continu- 
ation software package AUTO [13]. The bifurcation studies with time delay 
have been carried out using DDE-BIFTOOL [14]. 
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Fig. 1. Bifurcation diagrams as a function of the feedback strength k with and 
without time delay r. 



3 The Bifurcations 



We first carry out a detailed bifurcation analysis for our model system in 
terms of the feedback strength parameter k for different values of the time 
delay parameter r. Our results are shown in Fig. 1 where plots of the L2-norm 
(\z\) vs k are shown for b = —0.5 and r = 0.0, 0.5 and 0.6 respectively. The 
corresponding frequencies of the limit cycles are shown in Fig. 2 and one clearly 
notices frequencies tending to zero as the bifurcation point is approached. 
Going back to Fig. 1, we find that in the absence of time delay (r = 0) and 
for k > 0.42506 there exist two fixed points - a stable node and a saddle. 
At k = k c = 0.42506 we have an SNLC bifurcation such that for k < k c 
one only has a stable limit cycle branch whose period tends to infinity at 
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Fig. 3. Typical trajectories for various values of r (r = : solid; r = 0.5 : dash-dot; 
r = 0.57 : dash) illustrating a saddle separatrix loop bifurcation beyond a threshold 
Tt- The nullclines (N x ,N y ) along with the saddle (S ) and the node (N ) points for 
r = are also marked. 

the bifurcation point (see Fig. 2). The location of k c depends on the value 
of the shear parameter b but the nature of the bifurcation diagram remains 
the same as long as —1 < b < —0.5. The presence of finite time delay in 
the feedback introduces profound changes in the bifurcation diagram of the 
system. As seen from Fig. 1 the stable limit cycle branch now extends beyond 
the SNLC bifurcation point of k — k c and gives rise to a bistable region over 
this extension. The point k = k c now marks a saddle node (SN) bifurcation 



4 



point. The range of extension increases as a function of r. The existence of the 
limit cycle in the region beyond k = k c has a threshold character in r i.e. for 
a given k beyond k c the limit cycle appears only when r exceeds a threshold 
value (r t = 0.3759). Some typical trajectories for different values of r and a 
fixed value of k > k c are shown in Fig. 3. For r values below a threshold all 
trajectories end up on the stable node and there is no oscillatory behavior. 
At a certain critical value of r one notices the emergence of a separatrix loop 
trajectory and the consequent onset of periodic behavior. The dynamical origin 
of this behavior can be traced to the onset of an SSL bifurcation [1] occurring 
as the parameter r is varied. In both the cases of bifurcation (SNLC and SSL), 
the period of the limit cycle becomes infinite or the frequency becomes zero 
at the bifurcation point as we have already seen in Fig. 2 for the cases shown 
in Fig. 1. These bifurcations with their scaling laws have been beautifully 
explained by Strogatz [15]. 

The two parameter bifurcation diagram in k and r space is shown in Fig. 4. 
Different bifurcation branches as well as the regions of stable limit cycle/fixed 
point /bistable solutions are clearly demarcated. A system close to the meeting 
point of these three regions can exhibit four diffrent type of bursting [1]. A 
system can switch over from a steady state to the oscillatory state either 
through a saddle-node on a limit cycle (SNLC) or a through saddle-node (SN) 
bifurcation whereas the system can switch back from the oscillatory to the 
steady state through SNLC or a saddle separatrix loop (SSL) bifurcation. In 
principle, this gives rise to the possibility of four different types of bursts 
to arise: SNLC/SNLC (Type II or parabolic bursting), SNLC/SSL (no name 
yet), SN/SNLC (triangular bursting), and SN/SSL (Type I or square-wave 
bursting); the first bifurcation refers to the steady state to the oscillatory state 
and the second one refers to the oscillatory to the steady state bifurcation. 
The parabolic (SNLC/SNLC) and the square- wave (SN/SSL) bursting are 
quite common in neuronal modelling. We now check the feasibility of bursting 
behaviour of this model through a simulation work as described in the next 
section. 



4 Simulation of bursting 

For the purpose of simulating bursts, we add a slow stimulus to the model 
(Equation 2) which can drive the system through bifurcations of equilibria 
and limit cycles. The stimulus chosen is s(t) = ee int where e is the amplitude 
and Q is the frequency of a periodic signal. The amplitude and the frequency 
of s(t) determines whether the system evokes no-spikes (subthreshold signal), 
single spikes, or bursts. We choose e = 0.02 and ft = 0.01 such that the system 
gets across the bifurcation regions for all the four different feedback strengths 
k chosen for the present simulation. The other parameter values are b = —0.5 
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Fig. 4. Stability diagram in the parameter space of k and r. Note the various bifur- 
cation boundaries demarcating the different stability regions. 

and uj — 1. The time delay r = 0, 0.3 and 0.5 for Figs. 5,6, and 7 respectiveley. 
The four feedback strengths (k's) chosen are such that k = 0.41 and 0.42 are 
in the oscillatory region and k = 0.43 and 0.44 are in the fixed point or the 
rest region. The stimulus added plays the role of driving the system across the 
bifurcation boundaries. 



4-1 Bursting with no delay 



We first simulate the case with zero delay (r = 0) and the results are shown in 
Fig. 5. The bottom curve correspond to the case when k = 0.44 which is quite 
far from the SNLC boundary or the critical k c and the stimulus just manages 
to generate a single spike burst. As k is decreased to a value of 0.43 we get 
the parabolic bursts. Since the onset as well as the termination of the burst 
is due to SNLC, the period gets larger on both ends of the burst. The larger 
spacings on both ends of the burst is an indicator of that. The top two curves 
for k = 0.41 and 0.42 correspond to the case when the system is already in 
the oscillatory region and the stimulus pushes the system back and forth into 
the fixed point region to generate the bursts. The two cases, the one in which 
the system is in the fixed point region and the other in which the system is 
in oscillatory region to start with, are distinguishable by the fact that the 
system has shorter duration bursts in the first case and larger duration bursts 
in the second. In the first case the bursting activity correspond to a continuous 
state of quiescence or fixed point interruptedd by a spike or a cluster of spikes 
whereas in the second case it correspond to a continuous nonuniformly spaced 
spike train interrupted by short spells of quiescencee. 
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Fig. 5. "SNLC/SNLC" (parabolic) bursting in the absence of delay (r = 0.0) 
for four different values of the feedback strenght k. Other parameters are 
b = -0.5, u = 1, e = 0.02, Q, = 0.01. 

4-2 Bursting with finite delay 

We now delay the feedback so as to be still in the monostable region (i.e. 
r < r t ), and keep other parameters the same as in Fig. 5, the results are 
depicted in Fig. 6. One clearly notices an enhanced spike rate and an increase 
in the number of spikes per burst. These features get further enahanced if we 
increase the delay to get into the bistable region. Fig. 7 is for r = 0.5 with 
other parameters the same as in Fig. 5. The onset of the spiking is on the 
saddle-node off the limit cycle and hence spikes with a finite frequecies not 
tending to zero (see Fig. 2) whereas the termination of the oscillatory state 
is through an SSL or a homoclinic bifurcation which slows down the spike 
rate on this end of the burst. One does notices these features in the simulated 
bursts. An increase or a decrease in the time delay within the bistable region 
too correspondingly enhances or reduces the spike rate within a burst and its 
duration. 



5 Conclusions 

We have proposed a mathematical model for the spiking-bursting activity of 
a neuron. Eventhough the model is not based on physiology, it does simulate 
a number of features observed in neuronal bursting. When we examine the 
bifurcation scenario in two parameter space of the feedback strength and the 
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Fig. 6. "SNLC/SNLC" (parabolic) bursting with finite time delay (r = 0.3) but less 
than the threshold value Tt for the bistable region to appear. Other parameters are 
as in Fig. 5 
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Fig. 7. "SN/SSL" (square- wave) bursting in the presence of a large enough time 
delay (r = 0.5) to be in the bistable regime. Other parameters are as in Fig. 5 

time delay, we find the relevant bifurcations responsible for different types of 
bursting activities. The dynamics of our system is capable of generating highly 
nonlinear stable limit cycles, which mimic the spiking/bursting activity of a 
neuron, and a stable fixed point, which corresponds to a rest state. The sys- 
tem, more interestingly, also has a parameter regime where both these stable 
regimes (i.e. the oscillatory and the rest) coexist. Existence of such bistable 
regime is the reason for generation of square- wave bursting through hysteresis. 



The monostable rest and the oscillatory states, separated by a saddle-node bi- 
furcation on a limit cycle, enables the generation of the parabolic bursts. The 
spike rate within a burst and its duration (or the width) gets enhanced by 
a time delay which is below threshold so as not to push the system in the 
bistable region. A large enough time delay pushes the system into the bistable 
region and we obtain square-wave bursts with enhanced spike rate as well as 
the width of the square-wave bursts. 
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